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ABSTRACT 


Interaction between an electron beam and a fully ionized plasma has 
been studied with a view towards its application ina structure-less 
traveling wave tube. Three basic approaches, of varying degree of rigor, 
to the problem have been pursued and analytical solutions for the circularly 
symmetric case obtained. Comparisons between the methods of analysis 
are made. 

The writer wishes to express his appreciation for the assistance and 
encouragement given him by Professor GlennA. Gray of the U. S. Naval 


Postgraduate School in this investigation. 
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1. Introduction 

Investigations into wave propagation in plasmas have been conducted 
by a number of workers (1), (2), (3) with a view towards employment in 
plasma diagnostics and traveling wave tubes. Extensive bibliographies are 
contained in these works and the reader is referred to these for previous 
work in the field. 

Basically, the problem is a boundary value problem in cylindrica] co~ 
ordinates (r,98,z) as depicted in Fig. 1 witha static magnetic field parallel! 
to the z-axis. Throughout this paper, the case of coupling toa helix will 
be considered as this is Aa common method of coupling into or out ofa 
traveling wave tube. The mathematical model used for the helix is the 
simplest of those that have been developed, that of the "sheath helix". 
The reader is referred to any standard work, such as Pierce (4) for details 
on this point. The results will be used without comment in this paper. 

The mathematical model used for the plasma will be that customarily 
used in studies of this nature unless specifically stated otherwise. This is 
the "cold" plasma wherein effects of collisions, recombinations, neutrals, 
and thermal motion are ignored and the electrons are considered as forming 
a "cloud" against a background of positive ions which provide overall (DC) 
neutrality of charge, but do not otherwise appreciably contribute to the 
problem due to the relatively large mass of the ions which renders them 
practically stationary. 


Three approaches to the problem are presented and compared. These 
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are presented in order of increasing complexity, if not rigor, and the results 
are compared using the third method as the standard of comparison 

The first method is an extension of the works of Boyd, Gould and 
Trivelpiece and has, as the major distinguishing features, the assumptions 
that the magnitude of the longitudinal propagation constant is much greater 
than that of free space and that the electric field may be represented as 
derived froma scalar potential. This is Trivelpiece's “slow wave" or 
"quasi-static" approximation. This method has the advantage of simplicity 
albeit at the price of rigor, but the extreme simplicity alone is of considerable 
practical value provided, of course, that the results provide a reasonable 
approximation to the true case. A feature of this "slow wave" approximation 
is that a TE mode of propagation is denied by the first assumptions. While 
of no great consequence when dealing with drift tubes or waveguides, this 
is troublesome when one attempts to derive expressions for interactions with 
a helix. Matching boundary conditions at the "sheath helix" requires both 
TE andtTM types of solutions as is stated in Hutter (5) and can quickly 
be demonstrated. Thus, a dilemma presents itself. In attempting to extend 
this simple method, a free space TE solution will be assumed within the 
electron beam and plasma regions and the inconsistency ignored. 

The second method is an extension of the works of Rigrod and Lewis (6), 
and Brewer (7), the latter being essentially a generalization of the former. 
This method, as developed for electron beam studies, solves Maxwell's 


equations in a region containing charge, and, through a perturbational 





approach, takes the effect of the charge into consideration in the boundary 
value problem by replacing the rippled beam by an equivalent smooth beam 
with a surface current density. Brewer's model for the electron beam is 
simply a beam of electrons and does not postulate positive ion neutrali~ 
zation of the beam as some other common analyses do. The plasma, in 
this type of analysis, is treated as the limiting case of an electron beam 
ina plasma with zero charge density. Brewer, in his paper on the subject 
(7), obtains only a TM solution although a TE solution is not negated as 
in Trivelpiece's analysis. In this paper, an approximate TE solution is 
obtained, and interestingly, is shown to be coupledto the TM solution 

in such a manner that, if the TE solution is identically zero, then the TM 
solution is also zero, with the converse not necessarily true. 

The third method of analysis is considered the most rigorous and is the 
most complex, mathematically. It consists of solving Maxwell's equations 
in an anisotropic media using Kales' (8) method of solution. This method 
provides an exact solution of the mathematical model as described above 
and further amplified in appendix A. Interesting features of this method are 
the requirement of coupled TM and TE modes of propagation which cannot 
be zero independently, and a mode degeneracy with non-orthogonal modes. 

In this paper, the subscript o will be used to denote dc (static) 
quantities and the subscript 1, time varying quantities. Unless noted to 
the contrary, all quantities will be assumed to vary as 


% 





a = fur ar 
F = fo (49,3) + F-4%, DL whe 
with a z-dependence such that 


vi ~¥ F (9%) (1.2) 


Following the usual procedure, the factor e 


jut 


will be understood and not 


written explicitly. 


2. Method I, Trivelpiece's "Slow Wave" Approximation 


Following Trivelpiece (1), it is assumed that 
VXE, =-¢¥E, ~ O (2.1) 


which then allows us to represent the electric field vector as derived from 


a scalar potential, i.e. 
—E=--V? (2.2) 


From Maxwell's equations, we have 


ViDaV(EE)=VE-VO=0 2.3) 


where € is the tensor dielectric constant given by 
é;) 4 Ero (@, 
= = €!-#5a  &, O (2. 4) 
O O E231 


A simple derivation of equation (2.4) is given as Appendix A although it is 
worth noting, in passing, that the same result for the case of a plasma 
alone may be obtained from the Boltzmann transport equation (10) with far 
fewer restrictions upon the derivation. 


Proceeding formally, one then obtains 
ioe Ln é,,2 ene Eno +€ is 
= _ | rb +4 19, 28) + basa? NF 


(2.5) 


+ $518 se] =O) 





If a product type of solution is assumed, 


er 
b=ROO@L* 0.6 


the following differential equation is then obtained 


x6 4 124, +124, + fy *) = - (2.7) 


re rye nh 264% Ey) 


the solution of which is 


} = A, (T™) + EN, (Tr), |< cou(m 8) +Dainlnd) Lh Y (2.8) 


where A, B, C, and D are arbitrary constants. Following the usual pro- 


cedure in waveguide type propagation problems. this will be expressed as 


b = A Jon (Tx) +BY. (Te) pte es 


where 
T=? “gt (2.10) 
We then have 
E,, =~ [ATH Ca) + BTMy (Ta) | een ae 
afr? ty 
E, = Agent (Tr) + By jaMy (Ta) | 2 (2.12) 





{ne-tye 
A Yq, (Ta) + YM p(T )E" a 


where the prime indicates differentiation with respect to the total argument 
While it is implicit in the assumption behind equation (2.2) that the 
time varying magnetic field is essentially zero with respect to the eleciric 


field, the curl H equation 
VX Ht, = 7 wert, (2.14) 


will be used to obtain approximate values of the magnetic field componen‘s 
He will be set equal to zero, otherwise the existence of another mode of 
propagation would be allowed which would contradict the basic assumption 


equation (2.2). Taking the components of equation (2.14), one obtains 


QV 


J, IMe — 4 We, Es + 4 éabe | (2.15) 


=e 


tL 
JL OB 7 


an aH — (2.16) 
a <is 4we, | 4énF,, ¢ ene | 2.16) 


£5 (*4e)- 243 = Wey Ea3E yg (2647) 


We then obtain (where the primes again denote differentiation with respect 


to the total argument) 





EME, 4 Ej (AT In, (TR) +B TN (Ta)) = (718) 


oF 4M (AT (Tr) + BM,(T-))] 


Bp = 4p Gi (ATT) + BTM Ta) + 


Ea 4M [AJ (To + 8Mn(T)) (2.19) 
7 


From this point on, only the axially symmetric case (n=0) will be considered. 
Also, only the time varying field quantities are involved so the subscript } 


will be omitted to simplify the notation. 


E> = AY J, (Tr) + BYN (Tr) (2.20) 
ie = AT, (Tr) r BTN, (Tr) (2.21) 
Hy = dae 6, (AT To) r BTM (Ta)) (2.22) 





Be ~ ae ~46 (ATS (Ts) +B T™M,(T>)) 120233 


Equations (20) through (23) are the field quantities derived from equations 1] 
and (14), to be used in the boundary matching problem. 

As was mentioned in the introduction, it is impossible to matcha 7M 
solution alone to the "sheath helix" using the standard boundary condi- 
tions as given by Pierce (4) or Beck (9). Some form of the TE solution must 
also be used. For this analysis, a free space TE solution will be assumed 
with more to be said upon this assumption later. it is to be noted that. up 
to this point, no assumptions, other than those made by Trivelpiece (1) 


have been made. Fora detailed justification of these assumptions, the 


reader is referred to Trivelpiece's work. Extracting, the appropriate TE 
solution from Beck (9), the boundary value problem then becomes (referring 
Eo Fig. 1) 
RegionI OQSrs¢a (Electron beam within a plasma) 
From Appendix A, 
€, ¢ Eo 0 
E=G& 7% § O — 


Q O é; | 
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F224, tI, (T) li, = AgTo(pr) 


2 2S 
E Se ee 

Hy = ASE AT (Ti) Ea = EP AT, (p2) 
Region I] aSrsb (Plasma region) 

From Appendix A 

©, 4¢%g oO 
Ee = —~zé 
= Eo Gf 8 =I 0 [pas 


O O C33 

E, = B XG (Tan) + aM, (Tan) 

He= tert [8 Tad (Tar) + BaTaM, (Ta) 
Hz = Ag Tp (PR) 

Eg = LP AgT, (pr) 


Region III b*<r sc (Free space) 


E,= C,T, (Pr) + (Koy (pr) 

to = 40 C1, (pr) — a Ca ky (pr) 

Hy = C3Ty (pay + CyK (pa) (2.28) 
Eo = FL C3L, (pr) tat Cy Ki, (pr) 


i 





Region IV r2c_ (Free space) 


E, = Kk, (px) 


He = Ra D, K, (pst) 


(2.29) 
ty = Dako (pr) 
Eo = HM Ky (pr) 
The boundary conditions at the helix are: 
ae TE 
Eo = EE 
Eye — ~Es* cory (2236) 


Ht +Hé cory =HE + HE cory 


As the general problem is rather complex, several simpler cases will 
be considered first. These are: 
Case I, the case of an electron beam completely filling the interior ot the 
Plasma and the helix, i.e., a region J and IV problem. 
Case II, the case of an electron beam of radius a, less than the diameter 
of the helix, with the plasma filling the helix, i.e., a region [, II, and 
IV problem. 
Case III, the case of an electron beam of radius b, passing through a plasma 
of radius b, surrounded by a free space region and a helix at radius c, 
with free space outside the helix, i.e., a region I, [iI and l1V problem. 


Case IV, the case of an electron beam of radius a passing through a plasma 


m2 





of radius b. surrounded by free space and a helix at radius c with free 
space outside the helix, i.e., a region I, II, If], and IV problem. 

Determinantal relationships are obtained for all four of these cases 
with the algebra relegated to Appendix B- The determinanta! relationships 
are given below. 


Case I 
(2.31) 


ae + COURT (re) Kted 4 A GG IEC} — 9 
AVG Re (ed Yo (77 C) 
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eA" 
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(re)? 


(70) 


op 2 OM 
L (xd) ly 


(74) % 
pr 


TH (9) °L 


(a0) °y 7 a 
Ae10 2 7 


024)'9 (20) 'T 


(F4) Lf) (206)! (20)!'T 
(0%) °y 


+ (rN £9."9 braces : 


T (7%) 'C ERS "3 4 


Az 
Vieloo 


+ 


Orr 


W@W 3syj 


(IY (9%) IF 
(781) NF - 


(4°) 'N6%'TF 
(FE) 4 - 
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14 
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(| 
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-(0)'(ott) = H 


& 'S 
£E5 M5 








(011) °C fol) 
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3. Method II, (Extension Of Brewer's Method) 

This method is an extension of the method of Rigrod and Lewis .6) and 
Brewer (7). Attention is also invited to Beck's (9) excellent treatment from 
which this presentation proceeds. 

Let us consider an electron beam with no angular variation in charge 
density, fields or electron motion, i.e., 

4 


ers? 


We may then write the Lorentz force equation in cylindricai coordinates as 


i-né*=-% [Ex + Bg 76 | 3.1) 


£4 (nd6) =-Y ~4A + 5.3 (3, 2) 


f= ~ 7 [Fy ~ Bn? | oa 


provided we assume as does Beck and Brewer, that the contribution of 
Eig to the right side of (3.2) is small with respect to Bot. More will be 
said of this approximation later. Using Busch’s theorem from electron 


optics, we may write 


goo = 4. ES 7 Xe 6, (3.4) 


Where Ce is the initial radius of each electron and By the z-directed 
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magnetic field at that point (assumed constant). This may be stated as 


a 
ees Pen Je 


Substituting this into (3.1), one obtains 


v6 AL a 4 a 
Je =~ Ex ty pat — ut (3.6) 


Let us now perturb (3.5) by letting 


G= 0,1 FL = Ng T7y 


to obtain 


R 
= ae Ar 2° (3.7) 


Next, assuming that the DC ripple or scallop on the beam is small, let 


i FL St a a 
re =-W5, +3 hae 7 | 40 (3.8) 


Perturbing this expression as before yields 
Gps = za 3.9) 
(w+ 4 XU,)%— GP" | 
Perturbing (3.2) and disregarding the effects of AC magnetic fields on 


electron motion compared with electric field effects 


= i) E, , 
a (w+474, )a (3.10) 


Applying the continuity equation 


/, 





Q 


Wane ee gi 


Q (3.1]) 
to the AC charge density, one obtains 
- A40-F J 
tay e an n) +5 (3.12) 


(w+ged,)  —~ (w+pry,) = 2 


At this point, let us examine what has been done. The equilibrium 
relation (3.8) has been perturbed to obtain (3.9). As stated in the intro- 
duction, electrons are assumed to interact only through the electric field 
and the above relations are derived as if the individual electrons compose 
a stream with continuity of charge in this stream maintained through 
equation (3.11). If a second group of charged particles were present, such 
that the DC equilibrium condition remained valid, comparable expressions 
for the time varying quantities could also be written for the second group 
of charged particles. The model of the plasma set forth in the introduction 
fits these conditions quite well since overall DC neutrality from the plasma 
ions and electrons is maintained. Assuming that the continuity equation 
holds (which it must as no mechanism for loss of charged particles has 
been assumed in the model), one can write expressions comparable to 
(3.9), (3.10), and (3.11) (with U=9) for both ions and electrons, but, 
since the mass of the ions is so much greater than that of the electrons, 
their effects are extremely small at any frequency considerably above one 


megacycle per second and will be ignored. Designating the beam electrons 
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by subscript b and the plasma electrons by subscript a and dropping the 


subscript 1 from field quantities since only ac fields are to be treated, the 


following relations are obtained . 


4h 24; 
li = [ap bo a7 a) T e734 


(we gh, )a 


pigs flundiegen 
AG (Wt al Uy )*— Sez 


19 


,13) 


» 16) 


- 18) 


- 19) 





An. = 4H En 
Pea (y%— fee 3.20 


eo fen ng tT AL 


7L 
and, from Maxwell's equations in a media containing charge 
AU Pron ao Lae (3. 


From (3.13), (3.14), (3.19) and (3.20) 


_ 4 , 4% )y) 
Sn ee + BeHur HJ En. (3.23) 


( U)+4.8Up ) a 623 


AS 
NO 


It will be convenient to write 


(3.24) 
mom (3,22) and (3.23) 
—~ © df 4WMFE 
> pr on 2 ES (3.25) 


At this point, Beck shows that the second term of the last relationship is 
small with respect to the first for electron beams and may be ignored. 
This approximation is not as clearly well taken in this analysis and will 
not be made. Comment on the effect of making it will be made later. 
Solving for E. 
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7 y dF, 
EL P°+ 4a oye (3,26: 


Small DC beam scalloping has already been assumed and since the 
assumption that the magnitude of the perturbation is small is inherent 1n 
the perturbational approach, ro r. Further, if we assume that the dc 
magnetic field strength is everywhere constant and the cathode is not 


shielded, B ¥ B~ and 
O z. 


SO UR 
Ne 
wv C (3.27) 


Making this approximation and writing 
Tg = aaa. + Ritige + fyb 5 


ene Obtains 


Up. = & toa he + "he 749 Uo 
2 Wee, " (Use, yd 


¥4 4 Yo an yey 15) 
(4) +494, ) WT prey wk A - 5 aL 


Writing this as 


Jy L irs s(4 26) eae” 


(3.2 


mK 
(O 


and 


Zi 





T= FER =—_4— df 


AL Q (3.31) 
(Pr +gwuh) or ae 
and from Maxwell's equations in a media containing charge 


a 
ahi ies) 2 De aan Se a eae. 


there results 


1 2 fndh) + PE, = 0 a 
where 


An PUguee +l) Prgeph) 


“(Pr gum lh (40 +p) — Oh 


(3,34) 


The solution of this equation is, 


E, = © b($) +CaM, (¥.2) (3.35) 


The notation at this point has become rather cumbersome. For clarity of 
presentation and ease of manipulation, some arbitrarily defined constants 
must be used. Although the use of the components of the dielectric tensors 
as derived in Appendix A, is objectionable in that an equivalent dielectric 
is not being used in this analysis, these quantities are familiar to workers 
in the field and their use will facilitate comparison of this method with the 
other two methods presented in this paper. Introducing these expressions 


and a quantity 


In 





7 a) ; F 
ho P+ gBAl (3.36) 
the significance of which will appear Jater, we have 


Ff = 4W6,(6 -/) + ¥U,(E,-G)E (3.37) 


L = 4H & (E3~/) 


(3.38) 
~ X¥U,6,{6-6é 
fey ee? Slims 0G; u) 3.39) 
/ 
9. WH GAC aad 
iin 0 &3 
A ~~ P (3.40) 


AVE, B+ pPydy (G~Ey) — YF 


At this point, a convenient check on the development exists. !f we 


let Pb go to zero and make all the substitutions indicated in (3.34). 









we obtain 
y = | -p ain wy t— tu, & 
MD (3.41) 





UW, 
| a Ws — iy & 


The second term in the second bracket in the numerator can be shown to 


be absent if, in equation (3.22), the term containing I, is ignored as 


Les 





negligible. If this term is omitted, the result is then identical to that 


obtained by Method I (Trivelpiece's "slow wave" approximation). it 3 


i 


obvious that this term is not always negiigible, even for very slow waves 
ya Z. ‘ £ ® 0 ’ w ‘ “ ¢ ( ry 
(|4"|>>k"). Similarly, if the term containing i, is ignored in (3,22) 


and if we let pres 0, the results can be shown to agree with Beck's (9) 


confined beam development. 


Thus far, the TM solution has been obtained with the principal! approxi- 


mation being that E 


16 is negligible with respect to Bor “BZ in (3. 2) 


This, in effect states that the coupling of the TE mode, from which Ee 
is obtained, to the TM mode is negligible. An approximate TE solution 
which does not ignore the coupling of the TM to the TE mode wil! now be 


derived. From Maxwell's equations in a region containing charge, we 


obtain 


to (-28) + (A ty") HL = — p25) (3.42) 


| 


nee |.) 
th = at, | th + J, | (3 43) 


(3,44) 


and, from Appendix 1 


24 





ae -~4.WE, — Ue, 
ie | Ws— we 


\] 


i a “¢ (W+g¥Yo) i — UW, a 


— UF 
Assuming 
Ue = 4a Yea, + oe et 
we obtain 
=? a plapgtdy) Wy Pd oe 
tu — Wu® {iieg¥a,) 
En [a — coor 


which, for brevity, will be written 


5, Ee + We & (mies 
Substituting (3.44), (3.26), (3.36), into (3.49), we obtain 


‘ -|8 = — (€,~ ) 26, 


2 


We now have 


ZS 


3. 48 


3.49 





(3.51) 
ee pW, Ep (6-!) ai $f nee 
f, 1. on =) 
where 
Ff, = Ai Jp (%, 7) a AN, (q.) (3.52) 


Equivalently, let us write 
5 A 4,/ 2 / 
pry . PHS — 2H — _¥phur€, (6-1 ; "£ | 
fe He poe tie — phi = VEG) | pl, | 
fp P, JL Pe JL 
where the primes denote differentiation with respect tor. 
A solution to this inhomogeneous Bessel equation may be obtained 


by letting 


Uy 


Hp = AYl(r) + Buy(n} (3.54) 


FL = CY, (n) 


(3.55) 


Making these substitutions in (3.53) 
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yY + 4 + YU, : pi 
oe AK + CVU, S (€-/)| 


= (3.56) 
A Y ‘ a. 
= Ay PUL O 
ARR+CY U6 (6 ll ie a 
For brevity, this will be written 


£64) se £43) ao (3.57) 


If we require that 


- . (3.58) 
(peerce ta _ 
i, 


as we must to be consistent with (3.55), f(U,) will then be identically 


zero. We may, without loss of generality, set B=C=1. We then have 
1 ue Q 
Ue +#—* — & Uy = O (3.59) 
ri 
the solution of which is 
Us = C, zie (2) mn ne (pr) (3.60) 


which then furnishes the reason for the somewhat peculiar choice of 
constants made in (3.36) and 
A=h =- Yo, (6) 
RR +X) 


(oly) 
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The complete solution is then 
He = GAG ON) + GBM (Hr) + GL(4N+GK ga) 0.62) 


As was Stated earlier, the plasma alone is considered as the limiting 
case of zero electron beam charge density. Hence we may write expressions 
for the plasma region alone which are identica! in form to the preceding 
with wh 70. The additional subscript p will be added to distinguish 
between the two regions. 

As stated in the introduction, a distinguishing feature of this method 
of analysis is the replacement of the rippled beam with an equivalent 
cylindrical beam with a surface current density. Expressions for the 


surface current density will now be developed. Following Beck's procedure, 


we write (at r=d) 


G 


ee 15 U4 (3. 63) 





Ge anon Ug tie 


< on te 4 (3.64) 


Using (3.15) and (3.16), we obtain 


G,= §4,(6 "Eo E (d) (3.65) 
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_ ade, 
a, = — (6-1) E (d) (3.66, 
From Maxwell's equations, we obtain 
: = 
Hy = = t+ (3. 67) 
_,wal) 
Eg = ta + | (3. 68) 


Using (3.50), (3.24), (3.26) and manipulating the various constants 


aie 
iS iar + — (3,69) 


w 


a eee + eee DE | 


A (3.2.7 0) 


Ho = Sli +z Ease) AM) a 
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a — [C3 Z, (B) - CK (pa) 4 


(3.72) 
Ww an 
i — y | AT) tAM (5,2) 


All the field quantities required for the boundary value problem have 
now been obtained. Now we shall establish the boundary conditions. 
From Stratton (11), we obtain the requirements that tangential E must be 


continuous and that 


—- 


AS X| Hy H, | = (3.73) 


which then provides the following set of boundary conditions 


Eo. = Fon By Aan 
(3.74) 
The boundary value problem is then: 
Region I (beam and Plasma) 
E, = A, Io (G%) Gi) 
he i — 
Hp = ~ 70M [+73 A, J, (4,2) (3.76) 
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ie =A, Fa Jo(%7 tA3L, (Rr) 


Ey = 4M 1h (DAT te) — 2H AL 


i, (8B 


Region Ii (Plasma) 


a 4,4, (Yn) + BNE) 

4 - — bp ge BIE 29) t BMA 

6 AWA 7 ii Ma) 

Hy = & epg t 2) + By boa) + Bd, Ra") + By KB) 


E,= gw 1d, Ub (6,-) om 
; it Ca 5 Hig”) $M, ) 


Ap 
- as (p ) = By Ky co) 


Region III (Free space) 


Fo = C, To (pr) + Cok, (pr) 


3 


(3.78) 


(3273) 


(3.80) 


(3.81) 


(3.82) 


(3.83) 





H, = a * [65 (px) —G aK ( “) (3. 84) 


H= Cat (pr) +GK, (pa) 


(3.85) 

ee 
Sa IG i (pr) — Cy K, (p>) (3.86) 

Region IV (Free Space) 

Fy cs D K, (pr) (3.87) 

Us 
i 0, K, (pr) (3.88) 
H, = DK, (pr) (3.89) 

— tw 

Ea = to Do Ky (pr) (3.90) 


As in Section two, the same four cases will be considered with the details 


contained in Appendix C. It must be noted that, since this method calls 


oi 





for a rippled beam, space for the ripple must be allowed. Cases one and 
two must then be considered as limiting situations where the free space 


region has shrunk to zero. 
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4. Method III, Solution of Field Equations by Kales Method. 

In Appendix A, tensor dielectric constants have been derived which 
take into account the effects of the electronic motion in an electron beam 
and plasma. Using the tensor dielectric constant, we may then solve 


Maxwell's equations 


VX F = ee 7A H (4.1) 


WC T+ 40d eS 


(4.3) 


oA 44) 


for a region with zero charge and current density, since the effects of 
electronic motion are already taken into account by the tensor dielectric 


constant. It will be convenient to change the notation slightly so that 


Gf, £62 0 
E= #5. & O DS SE (4,5) 


0 Oo & 
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the net effect of which is to include the factor €| in each of the Gy 


Kales (8) developed a method for the solution of Maxwell's equations 


in anisotropic media such that the dielectric constant was isotropic but 


the permeability, uw, was a tensor of the same formas (4.5). 


This 


method has been used by Stafford (13) to study resonance phenomena in 


a plasma column, and by Johnson (14) for the case of a plasma ina cylin- 


drical drift tube. In the presentation that follows, Stafford's 


will be generally followed. 


Proceeding, we may write 


VXE=YXE + QXCVE 
UXE= GXE + BEX 
VXH=UXH + &X(-1)H 


eh 7" “\ 
XH = YUXH Tiny 
Separating transverse and longitudinal components 
= A 
XEo= ¢ a 'ae. 
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notation 


(4.6) 


(4.7) 


(4.8) 


(4.9) 


(4.10) 





aX (-YE, fs = —4 yl, (4.11 
UXH = 4 UW &33 £5 ay (4.12) 


a XC-YH, -YH,) = 4 Ws, E, + We, OXE, “oe 


solving (4.11) for a, VR 
ZN th 4 WH 7 
and noting that, if F ae 


a, X (@, X Fe (a, FQ, - (2 a, )\F =-F (4.15) 


we may obtain 





= = 4uaA,z 
b= ~ EL aX —-\Y& (4.16) 


Substituting (4.14) and (4.15) into (4.13) and manipulating somewhat, 
we may obtain 


= 2 
AW Ee ( @XH.) = WH Eat. 44 Ev 2 KG, 


(4.17) 





+ 4h YE - YH, Ss YA, 





Taking the cross product of a. with both sides of (4.13) 


Ve 2 +V A, = 4W é,, (a, X E,) — Wes Ep (4. 


and substituting (4.14) and (4.15) into (4.18), we obtain 


+ whueal (4,xH,) = Ax [Br — eee | 





ae (4, 
aN 4. 
-gWen WE + 4WiM EH 
re ae 
At this point, it will be convenient to define 
2. a 9 
iS = tT WG, 4 
a UW YA 79. 
(4 
from which follows 
2 a x 


Using (4.17) and (4.19) through (4.22), we obtain 


CAH A, = =f XV - UH, + (Rive Kesey) (4. 


+G ty é9E» + Kv, 


a7 


520) 


,2)) 


23) 





Solving (4.11) for Hi. taking the cross product of a with both sides 


and substituting in (4.13), we obtain 


> — opens at 
oe WYER _AXYH = 4H Et Wey Qs, Xe (4.24) 
+ oe 


Substituting (4.11) into (4.18) and then using (4.24) and (4.20) through 


(4.22) we may obtain 


(KAY E, = Ve (-YK'E, + ArH) + 


(4.25) 


GZ AXY(wukH, EK E;) 


If the divergence of both sides of (4.25) is taken, the divergence of the 
second term on the right side of (4.25) can be shown to be zero by vector 


identities. We then have 


(Kee) VE, = ye (-¥K"E, + Ua Hy ) (4.26) 


Performing a similar operation on (4.23) yields 
Y 7,4 Te a _ ys 
(KA) VK, = ule WEE, ~K A, ) (4.27) 


We may write D= €: E in the following manner 


= eas. Be. | 
o> 5. ~ 4 Eg(22X&) + @ 636 28) 
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If we take the divergence of D and note that (from (4.10) ) 
“\ on 


we obtain 


writing 
a en gE, (4,31) 


and making the necessary substitutions results in 


—~  &9 E39 
Vi Be = — eo yale i z, o&z (4.32) 


Let us now write 


ViH=O= VA, Se (4.33) 


If we now substitute (4.32) and (4.33) in (4.26) and (4.27), we may, 


after some lengthy but routine algebra obtain 


Veh (KAS) H, a YWlern €23 5 =e (4.34) 
W 
2 eG 
Yea a ee tr Kn Sa. RF = O (4.35) 
! I 
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For brevity, let these be written as 


OH tal, t#E, =O (4.36, 


® 
VE Gaia oO (4.37) 


In order to obtain a solution to this pair of simultaneous equations, 
the artiface of assuming that both ES and H. may be expressed 4s 
linear combinations of two other functions of r and € shal! be used. 
The constraints between the constants that must exist to permit a solu- 
tion for the two new functions will then be determined. Having found 
these two functions, we may then solve for E. and He and use the 
uniqueness theorem to state that these are the solutions. 


We shall now let 


ES = 7Au | t+ BR Us (4.38) 


H = g,4 - Go, 49, 


(4.39) 


Substituting these two relations in (4.36) and (4.37) and manipulating, 


we obtain 
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(4,40) 


ae Ug Cetra) ad La (Cr to | | Ge 


AG Bo 


ys, T 4a| Bagel zs ards) 
fifa - Ai 


(4,41) 


+ OY, aqrba)-glerrdg | =n 
$.- £&, 


faewe require the coefficient of U. in (4.40) and the coefficient of U 


Z j 


in (4.41) to be zero, we will have two equations of the form 
VF + d°F = 10 (4.42) 
the solution of which may be written 
F= [T(e~) ae 4, ar)| _ (4.43) 


Setting these coefficients equal to zero and manipulating we may obtain 


" s cp +e, ee ig, +p 


| PR ‘j (4.44) 


4) 


ee 2 > —_ = 





Se = Pat dey — Batt _— 
a Py @a (4.45) 
Equations (4.38), (4.39), (4.44) and (4.45) comprise a system of six 
equations in eight unknowns leaving two relationships which we may 


specify arbitrarily. We could let p =P =1 » However the choice 


} 


Se gs ae 
Se, fy = 2 (4.46) 


will result in a more compact notation. Making this choice, we find that 


A’ —_ atr+crt ant 4 bed (4,47) 


where the upper sign is to be taken with subscript one. In terms of the 


dielectric tensor (4.5) 


ae Pa Q 
4 =a KV S +€33) ~KE, = G te, ~&3)—-k Ea) 
% i] 


Yo (4.48) 


at YA,» és) 


It is also found that 
2 a C) 
- 44 


6 7 | (4.49) 


Before proceeding to find the field quantities themselves, it is to be 


noted that the method breaks down when the expression in the denominators 
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in (4.40) and (4.41) is zero. It can readily be shown that this occurs when 


a 
(A- C) = —4bd (4.50) 


which yields the following requirement on % for this condition to exist 


2 


y= 7 = 63 ~ 9. \Fn S33 = “Fn 633 —(6,- B316( 6&3) &] (4.51) 


This condition will be discussed in Section five and for the present we 
shall only note that it occurs for certain particular combinations of the 
system parameters and that when it does we must examine the problem 
more closely. 


Using (4.38), (4.39), (4.46) and (4.49), we may write 


cs Ne 
ES — a A A*T (an) +A AM, (a)} ra (4.52) 


ha 4, 
c- 


G 
(4.53) 


ao ae : 
He=\ 4A 4A fa-c) Llan) tA « (4:-OM (»)) 
F ot Fe ee 
eel 

as the most general expressions for E. and He with the time and 2z 

dependence suppressed. We shall now limit our consideration to the 
axially symmetric case. 

Using (4.23) and (4.25), we may obtain the remaining field components. 


After some routine operations, we obtain 
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p— P aye 
(£54) I~ Phy ~ Wh («Ay (P priy+ (J-. 27),yHIN — wv cree Y | iG "4,7 


-V 


P _ «+ o 
ys ae V4 04) 


% r'6 
N3lmS + C30 — ORY (ie NY | ry 


(554) Le? 


3 5 4 _ : 
| . Many yt Yam iY = 0-7) 4 FLW (w2p'e is oy t = oH (.¥ Ly | 


faa 


¢ 


/ <7 (¥-,] 


a 74 


(ash) 5m 4 (5 oT (wp) WY + tgen,k oF irr 
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C 


O = 5} (7) (03- Fm /n— 


(% aH) 
(24). OD O°) L4.109 OMY i 





— 40%) 2e\ (> hye m - 


(7 - —p) 
(2'r) POD — (717) 2¢h109(20/) 





oe tet 


(lau WH) . Beam = 
A100) — Poe ar yt ir)? "C9 yey 17 i 


)-47), vn HA (2'r) le Ip ptor F+ (27) ea 


(0-H) artlra~_¥ | (OIE Ther F + (20) f 


roo) Nev (PIP A089 FF 
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We may now proceed to solve the boundary value problem as was done 
in Appendix B for Method I. Since the procedure is the same, only the 
results will be given and the details dispensed with. 

For case I, (beam and plasma filiing a helix of radius c, we obtain 
equation (4.58) as the determinanta! relationship. Before proceding to the 


remaining cases it is convenient to define 
2 
im A, — G 
= a | 
4. d (4.5 





Fl 
CO 
—— 


y KA wk” Ba, 
Ki A! | (4.60) 





MM; 


9 2 
VAMP — ROWE EK WE, 
Ki— #t (4.61) 





lL; = 


The additional subscript p will be attached to denote these quantities 
for the plasma region. Using these parameters, the determinantal relation- 
ships for cases II, III, and IV are given by equations (4.62), and (4.63) 


and (4.64), respectively. 
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5. Comparison of the Methods 

Trivelpiece's "slow wave" approximation is obviously the simplest of 
the three methods even though it is still so complicated that computer 
results are required for a clear insight into the system, a fact which de- 
creases the value of its comparative simplicity greatly. It has several 
shortcomings. First, it neglects the ac magnetic fieids in its basic 
assumption, which later has the effect of denying the possible existence 
of a TE type of solution. This leads to difficuities when attempting to 
match a field solution obtained from this method to any case where a TE 
solution is required, as ata helix. Equation (4.53) shows that a TE 
solution does exist and that it has poles, indicating that, regardless of 
how small the arbitrary constants, the fields associated with the TE mode 
will become appreciable in the vicinity of these poles; hence this approxi- 
mation cannot be used in these areas. One of these poles is a Zero of €4 
as given by equation (A.23) and indicates that the extension of Triveipiece's 
method fails in an area not predicted by itself. Equation (3.41) also 
brings outa significant difference with method I!. 

The extension of Brewer's method appears to improve the Situation 
somewhat at a great increase in complexity. This method, by neglecting 
Eg in equation (3.2), in effect, neglects the effect of the TE mode upon 
the TM mode but does predict a TE mode which is a function of the TM. 
The starting assumption in the analysis is that (neglecting the B, z term 


in (3.2)) 


20 





By” »? Ee q 


het 


Using the relations developed within Section 2 this requires that the ra‘io 


9. a | 
Bon Gel) | (Ww +30) — Ce | (5.2; 


Br (PPL?) W(t) ; 





be small. Examination of this relation shows that it has a second oarder 


pole at 


U) = -4 5 U, | hors) 


| os a 
without corresponding zeroes. It is to be noted that for /3 > & and 
the vaiues of Ud, and Up, usually found in traveling weve tubes the ratio 
(5.2) will remain small. Also, the B, z term in (3.2) has been negiected 


since B 


3, is very small. Looking at this term again shows that z hes a 


pole at the frequency given by (5.3). The situeétion is no! cleer and the 
determining factor is believed to be losses which have not been included 
in the analysis. Examining (5.2) further shows the possibilities for failure 
particularly for fast waves. However no simpie statement of these condi 
tions reveals itself. 

In general, this method is believed to be an improvement over Method | 
but, due primarily to the complex boundary conditions, it leaves littie. if 
anything to recommend itself in preference to method 21). 


Method III, the Kales’ solution has several interesting features. First 


a1 





coupled modes are predicted, i.e. , TE and TM modes are related to each 
other by a constant and neither can be zero independently at other than 
special combinations of the system parameters. Second, a two-fold mode 
degeneracy is predicted since we may take either of two values for the 
radial propagation constant. This aspect presents complications in the 
boundary value problem. 

Suppose that we try to match the boundary conditions using only one 
of the two possible modes. We would then find thet only one arbitrary 
constant would appear in the field equations for region I and matching 
boundary conditions would give us four equations in three unknowns. The 
problem would then be over specified and would aliow solutions for, at 
most, particular combinations of the system parameters. By taking both 
modes of each degenerate pair, we may match any physically realizeable 
boundary conditions. 

Very closely related to the above is the fact that the degenerate (de ) 
modes are not orthogonal. From Churchill (12), we may state the genera! 


requirements for orthogonality as 


a, U, ( a.) a A» 4,'(a) 


tf 
O 


5a (5.4) 
BU, (4) + 4y4/(4) = O 

where U isa solution of a Sturm~Liouville system. Since these are 
precisely the relations which exist between ES and Hy and Ho and 


Eo for each of the modes individually, it is obvious that the azimuthal 
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modes (M ) are orthogonal for the same S- Since each of the s, is 
associated with a separate Sturm~Liouville system the degenerate (s.) 


modes are not necessarily orthogonal. Evaluating an integra! of the form 


b- 
fo J (42) J (frjdr (5.5) 


will demonstrate that, except for unusual combinations of the constants, 
the modes are not orthogonal. 

In summary, it may be stated that the attempted extension of a rela= 
tively simple method (Trivelpiece's) develops inconsistencies with the 
exact solution of the model and, until numerical results can prove to the 
contrary, must be used with caution. The extension of Brewer's method 
appears to improve the situation somewhat but the complexity of the 
boundary conditions deprives the method of any advantage of simplicity 
and hence, has little to recommend itself in favor of the Kales’ solution. 

The Kales’ solution, while complicated in its derivation, is not more 
difficult to handle afterwards than any other. It has the distinct advantage 
of being an exact mathematical solution of the model, the only limiting 
assumption being the derivation of the dielectric tensor itseif. [tis 
readily demonstrated that the method fails at the poles and zeroes of the 
components of the dielectric tensor, and at a few other special combina= 
tions of the parameters such that the problem degenerates into a much 
simpler one. The failure at the singularities of the dielectric tensor is of 


no great improtance since the model also fails under these conditions, i.e., 
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infinite fields and propagation constants are not found in nature. This 
method also has the advantage that it offers the possibility of extenstion 
to include the effect of collisions between charged particles since the 
conductivity of a plasma is a tensor quantity of the same form as the 


dielectric tensor. 
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APPENDIX A 
DERIVATION OF TENSOR DIELECTRIC CONSTANT 

It is permissible to treat a plasma as described in the introduction as 
an equivalent charge-free dielectric whose characteristics vary as a 
function of frequency. Further, the equivalent dielectric constant fora 
plasma ina magnetic field is a tensor quantity because the electric field 
vector and the displacement vector are no longer related by a single multi- 
plicative constant. The elements of this tensor are calculated by adding 
the convection current density to the free space displacement current 
density and setting the sum equal to the displacement current of the 
equivalent charge free region. Two cases shall be treated, first, that of 
a plasma alone, and second, that of a plasma with an electron beam 
through passing through it. 
PIASMA 


Neglecting a Py, term as a second order quantity, one obtains 


= 
ee 


quent +pr = gwel A.) 


Using the equation of motion and neglecting effects of AC magnetic 


fields on electron motion 


die = 25 4 (APY) = -S(E AXE) yn 


the following are then obtained. 





2 (A.3 
‘WAI = S232 + UA, 
tO 6 my le FC (A.4) 
ee 
Wy Sas, ee 
$9 = rm Ay AS 


Solving for the components of velocity yields 


a se ins WE 1c ae 
=~) [bea + aa As 


ures 


ae wee Wee in 


=-7 & 
os (A.8) 
4 


Substituting these components of velocity in equation {A.1) and solving 


for the tensor, 


E,, # Erg 0 
E=6 |1-p5 ry O (A.9) 
O 0 & 


06 





Where 


_ Ue ap 
Oy = Ww, — Ww (A.11) 
a 
Q =) =e 
Se) un iy (A.12) 


PLASMA AND ELECTRON BEAM. 

Assuming that the interaction between beam electrons and plasma 
electrons takes place only through the electric field (in keeping with the 
model described in the introduction), one may solve the force equation, 
and equation (A.2) for the beam electrons and the plasma electrons separ- 
ately. Equations (A.6), (A.7), and (A.8) are the solutions for the plasma 
electrons. For the beam electrons, assuming all ac quantities vary as 


Wt-¢ Zz ; 
5, e the force equation may be written 


a = p(w gh) Ap = = -<lE + FXG | ie 


from which, the following are obtained 


| £ 
flwr eT), = oe Come (A. 14) 


a7 





4(W+47U.) Ye = -<F, + UU, 


iA ¥S 
(U/+49U,) WR = -—F 
A Bad a, yO Ip (A. 16) 
Solving for the components of the ac velocity yields 
Vr = 7) gla rbolE = He Ee (A177 
i Fj NO PP pelea) 
ge ( Wt 4¢U, )?— WU, (We, ) % We a 


eee aa (w+ YU.) Ce 
Vig = | eee — (anya zap 18) 


E 
i — 12 . 
‘3 ” 5 (cce, 7M) (A.19) 


The convection current density is now written as the sum of the current 
densities due to the beam and plasma electrons and one then obtains 
AWE, +BY + Ket = fvert 20) 
Solving for the components of the tensor 
E= 6 |l-4 &2 € O (A.21) 


/ 


O O €3 
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where 


1 / W,> - a* tue —(W+4VU, )4 (A 22) 


c 


a 
we [vp |, ee 
= we lye —w* |T Wr) |aF— waa, 2% 


=/- Wp _ hy 
- uy? (Ww +4xU, )* (A. 24) 
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APPENDIX B 
SOLUTION OF BOUNDARY VALUE PROBLEM, METHOD !?! 
CASE I (Region I and IV) 
Extracting the appropriate field components from Section two, we have 
Region I 
ee I H_ =A, 1, (pr) 


(B.1) 
= _., 1W a 
Hy ag 70 ee A, Bs I, (pr) 


Region IV 
E=B, K) (pr) Ho=B,K (pr) 
(B. 2) 


=—- IWE, =— |W 
Hy BL 5 © &, (pr) EB By : K, (pr) 


With the boundary conditions at the helix (r = c) 
Es = FY a 
ei = BE 
eae Conny ea 
He + HecoTY = He +HgcoTy 


Matching these boundary conditions 


A, 8G (70) = BK, (po) (B. 4) 


Ay 4, (pe) = by EF k, (pc) 


(B.S) 
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A, od, (T) =~ COTF mAs Ame (pe)| (B.6 





AyTo (pe) + COTY anes AT 7 (79 = 


Ba Ky (ee) + cory [~6, 4% Kj(pc) 8.7) 


Solving (B.4), (B.5), and (B.6) for A_, B,, and A, and substituting in 


ll 


(B.7); one obtains, after some manipulation 


1, corPvRLD (7) | Klee) _ pat 9 9, 
pc kK (pc) p* K, (FC) yA7 (Tc) — Do 


Making the approximation that 
o a 
prs-(Y+k) ¥-¥ px-¢s 


and using equation (2.10), the determinantal relationship 


a 
wel + corte pitiled [huled , 4461059 | = 0 is) 


PC E (pC Ko (p 5 In(h ©) 


is obtained. 
CASE II (Region I, II, and IV) 


The appropriate field components, from section two, are 
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Region | 








F, =A eG) W7e= Ay £, (p”) 
(B.10 
Gye, } 
Hy = + y AT, J (4) Fg = Ag aT, (pr) 
Region II 
E, = B43, (~) T 8s {My (Ty) 
UW & &y 
Hp = ¢ a (8,15. (2) + 8 T3.M,(T )| 
(B.11,; 
He, = Ag. Ty (p-) 
Ep = ~A, t2MT, (px) 
Region IV 
Fo = DK (pr) Hz = Dy Ko (pr) 
ep puller [ - ay (B12) 
@ = ~My FE, (pr) Fa = Oe Ki (pr) 


Matching tangential E and H components at the boundary between the 


beam and the plasma (r = a) 


A, YJ, (1,2) = B,8I5(%]2) + By XM, (G9) 


pare _ g¢wh ey + 
217, J, (Ta) = y [6,7 J (Ty) t Bo IM, (7,4) (B.14) 


and eliminating the AL coefficients, one obtains 
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B (ed  & Ta S(T ae 
11 JoCha) = =& 7 (74) 


B if (13.4) 6, la Ny 4) - 9 (B.15) 


213, (Ta) -& 7, I, (7/4) 


At the helix, boundary conditions (B.3) apply 


~A, #44 Tipe) = Dy EL ky (pc) 





(B.17) 
B, V5 (13.0) + 8M, (150) = Ay cory tL T, (pc) a 
Agto (pc) t coTY ae pelg, To Tile C) 87a MeO) = 
G 
MaKo (Po) — 0, Cord a “Ky (PY) (B.19) 


Solving the first three of these equations for D, ; Ay: and D, and substituting 


in the fourth; one obtains, after some manipulation 
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B, XI cI (ect) 4 Kyle YI(B) _ corVA"e 
ae LT, (Po kK — pt ha O) 


b 4° Keo) - ] , | mca Ta ed (pc) 
PTR Ge) RIOT P| + Ee | ORG aS + ee MTS 


Qiy ga a 
Es i Ex Ty M, (T9 ve Bkile9 piles yy cmcor9) = 0 (B.20) 


Again making the approximation that 
pwe-qt 


one obtains 


‘ rato y : ae Ky oo J] 
(B.21) 


{Mo (T) compat AG GM (G0)+ Kea = 
aE eT paket 7 {NG 33 MK Mt Aes N,(%30)}| = © 


Setting the determinant of coefficients of By and B., obtained from equations 
(B.15) and (B.21) equal to zero, the determinantal relationship is then 


obtained. 
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rs a (pc 
Ae, cane Cae righ 


MoT) _ af€u €33 V.(%@) Jo(Ta%) [un E32 Jj (3a) 
eo ~ foe heey! — E (Fa) Ve & Jha) , 


AGI. coPVAy Wee Kl) or ) = 
eo. ¢ al aM (hot BARS Wy(30) O (p.22) 


CASE III (Region I, III, and IV) 


The appropriate field components, extracted from Section two, are 


Region I (Beam and Plasma) 


Ex = ACS, 7") H, = Ag 4, (pa) 
WE» E see (B. 23) 
Hg = 1 ATG, (7,7) E, = ma ET (pr) 


Region II (Free space) 


eS SS ae ( pr) a6 ‘Gy Kp (pr) 


tg = does Cz, (pr) ~Cy K(pr)| 
(B.24) 
Hz = Gl (px) ci Cy isp (px) 
- ud} 2 
Ea = FETT, (pr) ~ GK (Pry 


Region IV (Free space) 


f= OY, K, (pr) Hz = Dy K,(p%) 





W Eo oe 
Hy = ~D, #5"? K, (pn) Ey = Dn ZAK (pa) (B.25) 
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Matching tangential E and H components at the boundary between the 


beam and free space at r=b 


Ay T,(p4) = GL ®) +G ko (pé) (B, 26) 


W 


(B27) 


From the two equations immediately preceding, the following is obtained 


-_ kK, (ps) 
—_ C. + C 0 f \ 
Ay 3 i ies (pé) (B. 28) 


























I, (ps) (B29) 
It is then concluded that C, is zero and that A, ~ C,- Matching E. 
and H, 
AV (G4) = GI, (p) + Qk (pe) (B30) 
| WE, & _ gwe _ 
dP“ AT I(T) = Ape [OM - Ok (es) a 
(B.31) 
and eliminating A, , the following relation is obtained 
4 [aed — SS Le, | hye) 4 Kes) |- 
JoCh4) = P&T 3, (7,4) (56 § p&na(na} (B.32) 
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At the helix (r = c), boundary conditions (B.3) apply 


CD, Cp) +O Kylee) =D, Ko (pc) (B.33) 


84h 657, (pe) = PEK (ps 








(B.34) 

CIolpo) + Ok, (pe) = FAG (peor (B. 35) 
(3To (pc) t COTP ApeG), Wg) = a (3K ip) = 

Do Kp (pe) + COTY 1% IX, (pe) a (B.36) 


Eliminating D, ; D, , and C, , the following relation remains (after consi- 


derable manipulation) 


a 2 
G, pee _ Arcor ‘| t| Keipe) | - 
Filpc) Ky (pg) 2” Ky (pe ‘ “NT (pc) & (oe) ; 


(B.37) 


Setting the determinant of coefficients of the Co obtained from equations 


(B.32) and (B.37) equal to zero, making the approximation that prw-3¢ 
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and using the relations between the Th and y , the following determinanta! 


relation is obtained. 


Ip (pc) | Tp #) . tL (ps) 
Ce) K (PD IL Jo O74) V6 J; (7 €) 


L, (ac) ACoTY|| Ke les) , £K, (pa) 0 
Tile) Ke) P*Ko (pe) Jo (G8) $5 J, "T-6)| 


(B.38) 


CASE IV (Region I, II, Ill, and IV) 
The field components are as given in Section two forall regions. At 
r=a, the boundary between beam and plasma, the problem is the same as 


for Case JI. The result is 





| Bias o (Ty) =e Ey €23 nay i 


Jo (7%) E, & J; (74) 
(B39) 
By | Me (Ts4) _ En€zz M(RA)| = 2 
BG) VSS JG4) 
Atr=b, the boundary between plasma and free space 
BY Jp (44) + & YN, (7; &) = Gi, (pb) tT Cs Ko (pd) (B.40) 
1% [BHA FMR] = 1G (pb) ~ Cod | te.41 


Substituting equation (B.39) in the above and defining 
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G>T Ta GA) P23 PT ATT 
oS 


(R42 
H = M,(%4) J, (ha) - z 42 NM, (Ta) Je (7) 2) B43) 
After some manipulation, the following is obtained 
q Fe (P5) &y pTR(W, (5.0) G ~ I(T 8)H) ~ YT, (8M (T8)6 - 
‘ruin + Gy IK. (pb) &y p mM (1 b)G “JBAH) + 
Pasfeesie Tesan)] = 0 
(B. 44) 


At the helix (r = c), the situation is the same as for Case ilI. We then 


have 
Tr, (p) — A*Co7 PK, (po) K, Ged) | ee 
«ls p™ K, (pe) |+ C | ke oa ° alle 


lf the determinant of coefficients of the a obtained form equations (B. 44) 
and (B.45) is set equal to zero and the approximation that px -g¢ 
is made and the definitions of the Th used, the determinanta! relationship 


may be expressed as 
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O = Hower sr (rx a Hoe - 5 (ru) (90 "9 rey 
(20) 'Z 7 (27)°Y ~o 7 (20) Ir 
(94°98) XX | oujey| | OYA ~ OPE 


g 


X {A (9%) 2 - sym) (+4)'y + (Hous. ~ (74) ) (yd) 2y Ng SE 9M F— 
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APPENDIX C 
SOLUTION OF BOUNDARY VALUE PROBLEMS, METHOD II 


Case I (Region I and IV) 





Extracting the appropriate field quantities from Section three 


Region I (Beam and Plasma, 0 < r S c) 


E, =A J, (HO (C.1) 
hh v3 

Hp = pay" Item AT (44) (C.2) 

H = ARTA tA; (ro) (C.3) 


es U/ ¥ 4 YE, (6-1) OY ~— 
Region IV (Free Space, r 2 c) 
b= 0k, Go (C.5) 
— {We 
He = 7D, Ki (po) (C..6) 


i 





= D., K, (pc) 


fae 


(C.7) 


i] 


Ey tn Dyk, (pc) 


(C .8) 
The boundary conditions are derived from the basic conditions (3.74) 
and the boundary requirements of the "Sheath helix" model. Brief commenti 
has already been made that this case must be considered as the limiting 
situation of what amounts to case three (to be considered) with the free 
space region shrinking to zero thickness. if we view the problem in this 


light, the boundary conditions are found to be 


- (C.9) 


va 


==. Cova 


a Ox (C.11) 


Forming the quantities required by the last two relations, using (3.26), 
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(3.65), (3.66) and the relations given above for region I, we obtain 


Hy a5 = A, es fos (GO) + AG (00 | FAI (4 ¢) (C.13) 


glad oe YY, A 2 
Hat Ge = A, ak al ME) Th rr) (C.14) 


Applying the boundary conditions given above, we obtain 


D, Ky Cpe) i py) (C.15) 


Aue ¥%, MH Fol) A Tit - Av SAAT (A) = FAK pe) (0 16) 


ACh +42) 





DK, po), 4h coT¥ K, (pc) (C.17) 


Ay acer UY Di XG) masts IC) | tA3Z, (pC) + 


A 
A aoe (SOA 6) VF (8) eory= 18 


Wy f? 
De K, (pe) — D #406 K (pc)coTY 


Using the same method as used in Appendix B, the determinantal relation- 


ship is found to be 
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gun th Weole NTO 4 EG T(g¢) + 
lh (As; 42) Gon 


44M L, (#,¢) Yok Ceghern Te) + ABS 
4 af, 


(C.19) 


no gaa im En) +B rH) THC) COTY + 
4 


EGIK 2) 4 sWleecoT¥EWAKC] = o 
UM wry K, (pe) PK; (pc) 


It is believed obvious, at this point, that the complexity of this 
method is considerably greater than that of Method I. It is also believed 
that nothing is to be gained by writing lengthy, complex expressions from 
which no insight into the nature of the problem can be obtained by in- 
spection. The method of presentation of the determinantal relationships 
will, therefore, be changed to expressing the relationship in the form of a 
determinant set equal to zero. 

Case II (Region I, II, and IV) 


Extracting the appropriate field quantities from Section three 
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Region I (Beamand Plasma, 0 Sr 


Ey 


Hz 


= 
ee, 


er Sp Cone) 


< a) 


Ay Fy Jol(Gy) tA3 Zo (27%) 


A de) FM AST (37) 


Region II (Plasma, 


E,= & 


Ho 


H, = 


= oe Pe 


6, Fin GL, 


asr s&s oc) 


nr) +8 FM 


44 


_ +% BT (mp) tBoM(¥ 
qu Re ON S72 awl et) 


7) tBLo (pn) +8y KC 
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20) 


(C21) 


(Gi; 


(Cc 


(C. 


© 


23) 





= At hp WeEo( En (40) tM (G7 
hy (hp t%) Gs * } 


“FL 8 (R) — Ey kh”) 


ap 


Region IV (Free Space, r 2c) 
Fs = D, i (pc) 
Hy = — 9, 42> K (po) 
ie = Do Gs (pe) 


Eg = Dae h, (po) 


(C.27) 


(C.28) 


(C.29) 


(C.30) 


(C.31) 


The boundary conditions for this case must be examined. At the 


helix (r = 0), the same considerations apply as were discussed in 


case I with the simplification that G, is zero. At the boundary be- 


tween the beam and the plasma, the situation is more complex. 


The 


boundary conditions expressed as equations (3.74) apply, but we must 


examine the expression for Go . We can write an expression simi- 


lar to (3.66) for the beam and also one for the rippled inner surface of 


the plasma region. Somewhat heuristically, the surface current 


arises due to the difference in the two media. Therefore, if we 
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let the difference between the two media approach zero ina smooth manner. 


this surface current must also approach zero. Gg must then be 


Ons ——— — eee ne (C.32) 
After some manipulation, we obtainat r=a 
Pere (C233) 
b, 6 icin FE ( ao n) th, yr J Un2) 
ha &p =| Fils (a) + Htayn(6-1 To) + 
A 
Ay Telpa) + B, RbpUAlENET (4 a) + 

ae (C,34) 


Bo i ATE N, (bn %) 
Be 


andat r=c 


Hebe = 8, re tt,9 7 Pe CELE) ede Tbe) + (6.35) 
Ap 


Bs pe } byt 6M + 
27 


Ap 


Bs T(z ¢) + By Ky (A, ¢) 


a 





By applying the appropriate boundary conditions at the interfaces 
between the regions, we may write eight equations in eight unknowns 
(the arbitrary constants) and write the determinantal equation as the con- 
dition for a non-trivial solution, i.e., the determinant of the matrix of 
coefficients must equal zero. When this is done, equation (C.36) is 
obtained. 

Case III (Region I, III, and IV) 


Region I (Beam and Plasma, O<r <a) 


EF, =A, Je (Mr) 


(C.37) 
He = a tbe Ay TJ) (br) (C.38) 
hh = A ATA) FAST (pr) (C.39) 
Ey = 


A4MIe, KO (6-) A T¥,r) - | 
BR +E) er eta") (40 


H2(a)-G, (a) = A, as (EY J(ta) +B Ie “a FA35,(82) (¢,.41) 
2 
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Ay) 0 “yy? ve 1% 
| hail ya Ody = O 7) Ay (0 7) ae - 1%4)'M(/-"9) 99° ER oe Ale (1-9) 05 Age QO QO 
Pe ee c Cie 
(0%) WOn"5)97me FA * OPAL (L'9)29 MRA 
tuntlly ae (ody 0 0 O O O 0 
“ (Tht) (i GY 
: 4% / Jf 3% / of A r7/% wie 
(Mp =O OO ae — OO T  loailestsindaenn® OI y nan Y J 
Q A. 
O (04) °y- O O (O°AW Cah O O 
Cree 
4 yyy 
ae | vp'y _¢ Atyd ) 19 a 
(4 7) ) 0) ( eae ve (-'3) me ent ie 
lee a a, ? 
oe 
| O O (nny 0% °r (0 ee ae 
Thay °c (i-'9)°a mp ut 
dj of 44(9- agree 
0 O 0 . v csv] he JE et | 4 | “ Nani 
| aol X ATE 


O O 0 0 (oA - (of - 0 fox) 





Lh we ; Qn, 2 
HA) +f) = 4, cn ¢ MA (6° €y) (7én) th te yl WA) (6.42) 


Region III (Free Space, a =T1r < c) 





FE, = CZ (pr) +, (pr) (C. 43) 
_ 4Wls 


Hp = C3 To Pr) + Gy ky (pr) 





(GaAs) 
._ —4w 
Ey ~ oe Gx (pr) ar, k (en (C. 46) 
Region IV (Free Space, r @ c) 

FL = D, Ky Cp») (C.47) 

__ pn {WE 
eer DF Ky (p>) (C. 48) 
H2 = Do Ko lpr) (C. 49) 
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_ Ue | 
Ep, = 7 1K) (pr) (C.50) 


Applying appropriate boundary conditions, as in Case II, the 
determinantal relationship is found to be equation (C.51). 
Case IV (Region I, II, III, and IV) 

Examining this case, we find that all the required quantities have 
been developed for the preceding cases, therefore equation (C.52) is given 


as the determinantal relationship without further comment. 
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